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SEMICLASSICAL AND SPECTRAL ANALYSIS OF OCEANIC WAVES 



CHRISTOPHE CHEVERRY, ISABELLE GALLAGHER, THIERRY PAUL, 
AND LAURE SAINT-RAYMOND 



Abstract. In this work we prove that the shallow water flow, subject to strong wind forcing 
and linearized around an adequate stationary profile, develops for large times closed trajec- 
tories due to the propagation of Rossby waves, while Poincare waves are shown to disperse. 
The methods used in this paper involve semi-classical analysis and dynamical systems for 
the study of Rossby waves, while some refined spectral analysis is required for the study of 
Poincare waves, due to the large time scale involved which is of diffractive type. 



< 

^ ■ 1. Introduction 

The problem we consider is motivated by large-scale oceanography: the main physical phe- 
nomenon leading this study is the existence of persistent oceanic eddies, which are coherent 

ff^ . structures of vortex type, spreading over dozens of kilometers and propagating slowly over 

^ I periods from one year to one decade. These structures have been observed long past by 

physicists [13 HZl [THl [2U1 EZ] who gave heuristic arguments (reproduced below) to explain 
their formation due both to wind forcing and to convection by a macroscopic zonal current. 
Giving a (much less precise) mathematical counterpart of those arguments, even at a linear 

l/^ ■ level, requires careful multiscale analysis and rather sophisticated tools of semiclassical and 

^^ I microlocal analysis. In this paper we simplify the model by considering particular macro- 

scopic currents, which are stationary solutions of the forced equations. This allows to exhibit 
trapped Rossby waves, by solving the dynamics associated with an appropriate integrable 
Hamiltonian system. We prove also that the other waves produced by the dynamics, namely 

^ ' Poincare waves, disperse on the same time scales (which turn out to be of diffractive type). 



1.1. Physical observations. Simple observations show that large-scale ocean dynamics 

can be decomposed as the sum of the solid-body rotation together with the Earth, convection 
by macroscopic currents (such as the Gulf Stream in the North Atlantic, the Kuroshio in the 
North Pacific, Equatorial or Circumpolar currents), and motion on smaller geographical zones, 
due for instance to the fluctuations of the wind and more generally to the coupling with the 
atmosphere. While the spatial extent of macroscopic currents is of the order of a hundred to 
a thousand kilometers, those fluctuations are typically on dozens of kilometers. We therefore 
expect eddies to be particular forms of those fluctuations. The point is to understand why 
they are quasi-stationary, or in other words why they do not disperse as other waves. At this 
stage we have to describe briefly the different kinds of waves that can be found in the 
ocean as linear responses to exterior forcing. They are usually classified into two families, 
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depending on their typical period and on their dynamical structure. The exact dispersion 
relation of all these waves can be computed explicitly (U [101 El El E] in simplified cases (no 
convection, linear approximation of the Coriolis parameter). 

• Poincare waves, the period of which is of the order of a day, are fast dispersive 
waves. They are due to the Coriolis force, that is to the rotation of the Earth ; 

• Rossby waves propagate much slower, since the departure from geostrophy (that 
is equilibrium between pressure and Coriolis force) is very small. They are actually 
related to the variations of the Coriolis parameter with latitude. In particular, they 
propagate only eastwards. 

The heuristic argument leading to the existence of quasi-stationary coherent structures is 
then as follows (as suggested by physicists): the wind forcing produces waves, in particular 
Rossby waves which would propagate, in the absence of convection, with a speed comparable 
to the bulk velocity of the fluid v ~ lOrns"^; the convection by zonal flow may then stop 
the propagation, creating ventilation zones which are not influenced by external dynamics, in 
particular by continental recirculation. We are then led to studying wave propagation under 
the coupled effects of the pressure, the Coriolis force and zonal convection, that is to studying 
a system of linear PDEs with non constant coefficients. 

1.2. The model. The system we will consider is actually a toy model insofar as many physical 
phenomena are neglected. Our aim here is only to get a qualitative mechanism to explain 
the trapping of Rossby waves. More precisely, we consider the ocean as an incompressible, 
inviscid fluid w^ith free surface submitted to gravitation and wind forcing, and further 
make the following classical assumptions : the density of the fluid is homogeneous p = po = 
constant ; the pressure law is given by the hydrostatic approximation p = poQZ ; the motion is 
essentially horizontal and does not depend on the vertical coordinate, leading to the so-called 
shallow water approximation. For the sake of simplicity, we shall not discuss the effects 
of the interaction with the boundaries, describing neither the vertical boundary layers, known 
as Ekman layers, nor the lateral boundary layers, known as Munk and Stommel layers. 
We consider a purely horizontal model, and assume an infinite domain for the longitude 
(omitting the stopping conditions on the continents) as well as for the latitude (this may 
be heuristically justified using the exponential decay of the equatorial waves to neglect the 
boundary). The evolution of the water height h and velocity v is then governed by the 
Saint- Venant equations with Coriolis force 

dt{poh) + V-{pohv)=0 
(1.1) 

dtipohv) + V • {pohv (g) f ) + uj{pohv) + poghVh = p^hr 

where ui denotes the vertical component of the Earth rotation vector 0, v := {—V2,vi), g is 
the gravity and r is the - stationary - forcing responsible for the macroscopic flow. It depends 
in particular on time averages of the wind forcing, temperature gradients and topography. The 
equations are written in cartesian coordinates (xi, 2:2), where xi corresponds to the longitude, 
and X2 to the latitude (both will be chosen in R). The vertical component of the Earth rotation 
is therefore r2sin(x2/i?), where R is the radius of the Earth, but it is classical in the physical 
literature to consider the linearization oi uj (known as the betaplane approximation) uj{x2) = 
QX2/R', ixLost of our results will actually hold for more general functions co, but in some 
situations we shall particularize the betaplane case in order to improve on the results. In 
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order to analyze the influence of the macroscopic convection on the trapping of Rossby 
waves, we will consider small fluctuations around the stationary solution 

h = h, V ■ {v ® v) + ujv = T, div v = 0, 

where h is a constant. Physical observations show that the nonlinear convection term is 
essentially negligible compared to the Coriolis term, so that the previous equation is nothing 
else than the Sverdrup relation (see |27j ) . 

1.3. Orders of magnitude and scaling. Let us introduce the observation length, time and 
velocity scales ^o (of the size of the radius of the Earth R), to and vq, and the nondimensional 
variables x = x/Iq, t = t/tQ, and u = {v — v)/vq . We also define the typical height variation 
5h and the corresponding dimensionless variable rj = {h — h)/6h. We denote by Vc the typical 
value of the velocity of the macroscopic current: u = (ni,0) = v/vc ■ The length scale /q and 
the convection velocity Vc are fixed by the macroscopic flow: typical values for the Gulf Stream 
are Iq ~ 10^ km and Vc ~ 10ms~^ . As we are interested in structures persisting during many 
months, a relevant choice for the observation time scale is to = 10® s (~ 0,38 months) . The 
associated Rossby number is then Ro := l/(to|r2|) = 0.01, recalling that |ri| = 7.3xl0~^s~"'^. 
The variations of water height which can be observed are typically of the order (5/i ~ 1 m to be 
compared to h ^ 10^ m . The influence of gravity (through hydrostatic pressure) is measured 
by the Froude number Fr := {vqIq) / {t^gSh) ~ 0.1, considering namely fluctuations of 
order vq ~ O.lms"^. Defining e := Fr and dropping the tildas (note that as often in Physics, 
e is not really a very small number), we therefore end up with the following scaled system 



(1.2) 



dtr] + -V ■u + u-Vrj + e^V • (r/u) = , 



e 



dtu H — ^-ftu"*" H — Vr/ + u ■ Vn + u ■ Vu + e^u ■ Vu = , 



£2" e 



where h := a;/|il|. We shall compute the response to the wind forcing, assuming that the wind 
induces a pulse at time t = 0: since the wind undergoes oscillations on small spatial scales, 
the initial data is further assumed to depend both on x and x/e. Typically 

( k- 
(1-3) {rie,Ue)\t=Q = {'nk{x),Uk{x))exp H — 

for some k ^ 1?. More generally we shall consider initial data which are microlocalized (in 
the sense of Appendix IB]) in some compact set of T*'R?. 

1.4. Local well-posedness. The local existence of a solution to the scaled Saint- Venant 
Coriolis system ()1.2p supplemented with initial data in the form (jl.Sp comes from the general 
theory of hyperbolic quasilinear symmetrizable systems. Defining the sound speed uq by 

77=[(l + e\o/2)'-l]^^ 
we indeed obtain that (|1.2p is equivalent to 

(1.4) £^dtU + A{x2,eD)U + £^Q{U) = Q, U = {uo,ui,U2) 
where A[x2,£Dx) is the linear propagator 

''eu ■ eV edi £82 

(1.5) A{x2,£D):=\ edi eu ■ eV -b{x2)+e^u[ 

ed2 b{x2) eu ■ eV 
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and Q{U) := Si{U)ediU + S2{U)ed2U with 

/ Ui i^o \ / U2 i^o 

(1.6) SiiU) := ino m and S2{U) := U2 

V ni / \luo U2 

Because of the specific form of the initial data, involving fast oscillations with respect to x, 
we introduce semi-classical Sobolev spaces 

HI = {U eL^ /\\U\\hi < +00} with WUfffs = Y^ \\{eV)^U\\l2 . 

\k\<s 

We shall also need in the following to define weighted semi-classical Sobolev spaces (in the 
spirit of [To] ) , adapted to the linear propagator as explained in Section [Tj 

(1.7) W,^ := {/ G l2(r2) / (1 - e'df)Hl - e^dl + h\x2))^f G ^^(R^)}. 

A classical result based on the Sobolev embedding (see Section [7] for related results) 

C 

||eVC/||Loo < — ||Vf||//s for any s > 1 , 
e ~ 

implies that (fO|) has a unique local solution U^ G L°°([0,re), iJ|+^). Note that the life 
span of Us depends a priori on e. One of the goals of this article is to show existence on 
an e-independent time interval. 

2. Main results and strategy of the proofs 

Most of this paper is concerned with the analysis of the solution to the linear equation 

(2.1) e^dtV ^A{x2,eD)V = ^, V = (v^,vx,V2) 

which is expected to dominate the dynamics since we consider small fluctuations. The de- 
scription of the linear dynamics is provided in Theorem [T] below. The comparison between 
linear and nonlinear solutions is postponed to the final section of the paper (see Theorem [2]). 

For technical reasons we shall restrict our attention in this paper to the case of a shear flow, 
in the sense that u{x) = ('Ui(a;2),0), where ui is a smooth, compactly supported function. We 
shall further assume for simplicity that the zeros of ui, in the interior of its support, are of 
order one. We shall also suppose throughout the paper that 5 is a smooth function with a 
symbol-like behaviour: 

(2.2) Va G N, 3C„, Vy G R, |6(")(y)| < C7„(l + b\y)), 

and we shall further assume that lim b {y) =00, and that b has at most a finite number of 

critical points (that is to say points where b' vanishes). Without such an assumption one could 
not construct Rossby waves. We shall also suppose that the initial data is microlocalized 
(see Appendix[BJ) in some compact set of T*R^ (which we shall identify to R^ in the following), 
denoted C and satisfying 

(2.3) cn {6=01 = 0. 
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Thanks to this assumption, which is propagated by the hnear flow, one can diagonahze the 
system into Rossby and Poincare modes. Finally in order to avoid pathological trapped Rossby 
trajectories we shall also require that 

(2.4) cns = 0, 

where S is a codimension 1 subset of R^ defined in Proposition 14. 6[ 

2.1. Statement of the main results. In this paragraph we shall state the two main theo- 
rems proved in this paper. The first result deals with the linear system (|2.ip . 



Theorem 1 (The linear case). There is a submanifold A o/R^, invariant under translations 
in the xi — direction, such that the following properties hold. 



Let Ugfi be e—microlocalized in a compact set C satisfying Assumptions li2.3\) - (2^^ . For any 
param,eter e > 0, denote by Vs the associated solution to i2.1\) . Then for all t > one can 
write Vs{t) as the sum of a "Rossby" vector field and a "Poincare" vector field: Vs{t) = 
Vj^{t) + V/'{t), satisfying the following properties: 

(1) There is a compact set K ofB? such that 

Vt>0, \\V,''{t)\\L2^K^^0{en 

if and only if the e-frequency set of V^^(O) intersects A. 

(2) Suppose that 6^ has only one non degenerate critical value (meaning that (6^)' only 
vanishes at one point, where (6^)" does not vanish). Then for any compact set Q, 
in R^, one has 

Vt>0, ||y/(t)||i2(f,) = 0(e°°). 

In particular supposing that 6^ has only one non degenerate critical value, then there is a 
compact set K of R^ such that 

Vt>0, \\Ve{t)\\L2^K)^0{en 

if and only if the e-frequency set ofV/^{0) intersects A. 

Remark 2.1. Actually A corresponds to the set of initial positions and frequencies in the 
phase space giving rise to trapped trajectories for the Rossby hamiltonian. This will be made 
more precise in Section [^ where we shall prove that under some additional (non restrictive) 
assumptions on u, A is of codimension one. In particular it will he shown that some of those 
trapped trajectories actually exhibit a singular behaviour in large times, in the sense that they 
converge in physical space towards a point, while the ^2 frequency goes to infinity. This could 
he interpreted like the creation of some sort of oceanic eddies. 

The result ^ is related to dispersive properties of the Poincare hamiltonian on diffractive 
type times (of the type 0{\/£^)), which requires some spectral analysis. Due to the assumption 
on b^ one can write a rather simple proof; more general conditions could he treated, but the 
Bohr Sommerfeld quantization would require to decompose the phase space into various zones 
according to the geometry of the level sets of the Hamiltonian, which is much more technical 
and beyond the scope of this article. Actually in [H] we propose a different approach, based 
on Mourre estimates, which allows to relax very much the assumptions on b^ and on u. 
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The final section of tliis paper is devoted to the proof of the following theorem, which states 
that the very weak coupling chosen in this paper implies that nonlinear dynamics are governed 
by the linear equation. We also consider more generally the following weakly nonlinear system 
(with the notation ()1.5p and (I1.6P ): 

(2.5) e^dtU + A{x,eD^)U + e^+''Si{U)ediU + e'^+iS2{U)ed2U = 0, 7?>0. 

The case r] = corresponds of course to the original system (jl.4p presented in the introduction. 

Theorem 2 (The nonlinear case). Let Usfi be any initial data hounded in W^. Then the 
following results hold. 

(1) The case ij = 0: 

(a) There exists some T* > such that the initial value problem (2. 5\) with rj = has 
a unique solution Us on [0,T*[ for any e > 0. 

(b) Assume that the solution Vs to the linear equation Ii2. l\l satisfies 

lkK||L2([o,T*[;L«=) -^ as e ^ 0. 
Then the solution Us to h2. 5)) with r] = satisfies 

\\Us — V^||l2 — ^ uniformly on [0,T*[ as e — )• 0. 

(2) The case r] > 0: Le T > be fixed. Then there is eq > such that for any e < Eq, the 
equation Ii2.5\) has a unique solution Us on [0,T]. Moreover, 

\\Us — Vs\\l2 —7- uniformly on [0,T] as e — t- 0. 

Remark 2.2. Result J76)). joint with TheoremUl implies in particular that as soon as b^ has 
only one non degenerate critical value, then for positive times the energy of Us on any fixed 
compact subset is carried only by Rossby waves. The refined L°° estimate on the linear solution 
required in result J76]) should be proved by using WKB tools. For the sake of simplicity, we 
shall not consider such technical estimates here, all the less that we do not expect them to be 
enough to get an optimal result regarding the nonlinear problem (see Remark \7.1\ ). That is the 
reason why we consider, in result (0), a weaker coupling still. That result implies in particular 
that the L^ norm of Us on any fixed compact subset may remain bounded from below only if 
there are trapped Rossby waves, i.e. only if the e -frequency set of the initial data does intersect 
A (with the notation of Theorem\^. 

2.2. Some related studies. This work follows a long tradition of mathematical studies of 
fast rotating fluids, following [28] and [18]; we refer for instance to [^ and [TJ] for a number 
of references. The present study concerns the case when the penalization matrix does not 
have constant coefficients. A first study in this type of situation may be found in [12], where 
a rather general penalization matrix was considered. Due to the generality of the situation, 
explicit computations were ruled out and no study of waves was carried out. In order to 
compute explicitly the modes created by the penalization matrix, various authors (see [^i [Sj, 
[lOj as well as [13]) studied the betaplane approximation, in which the rotation vector depends 
linearly on the latitude. In that case explicit calculations may again be carried out (or some 
explicit commuting vector fields may be computed) and hence again one may derive envelope 
equations. In this paper we choose again to work with a more general rotation vector, this 
choice being made possible by a semi-classical setting (see the next paragraph); in particular 
that setting enables us diagonalize the system approximately, therefore to compute waves; 
note that a related study is performed by two of the authors in [5] via a purely geometric 
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optics approach, where no exphcit diagonahzation is performed (actually the initial data is 
strongly polarized so that only Rossby modes are present, including in the non linear setting) 

Another feature of our study is that it is a multi-scale problem, in the sense that the oscillation 
frequency is much bigger than the variation of the coefficients of the system. This is dealt 
with by using semi-classical analysis (which, compared to the previous paragraph, enables 
us to compute almost commuting vector fields although the penalization matrix no longer 
depends only linearly on the latitude). Such techniques are classical in geometrical optics, 
but in our case the additional difficulty is that the propagators are linked to different time 
scales: this is due to the fact that the system has eigenvalues at different scales (one is 
actually a subsymbol). In particular we are mostly interested in the role of the subsymbol 
in the dynamics, as this subsymbol is responsible for the trapping phenomenon we want 
to exhibit: this implies, by semi-classical analysis, the need to study the dynamical system 
induced by that subsymbol. On the other hand this also means that the dynamics linked to 
other eigenvalues must be analyzed on diffractive-type time scales, therefore much longer than 
that allowed by semi-classical analysis. We are able to show the dispersion of those waves by 
using spectral analysis and Bohr-Sommerfeld quantization. 

2.3. Organization of the paper. The proof of Theorems [1] and [2] requires a number of steps 
which are described in this paragraph. 

2.3.1. Reduction to scalar propagators. Persistent structures are related to the propagation 
of Rossby waves. Our first task is therefore to transform the original linear system ()2.ip into 
three scalar equations. One is polarized on Rossby waves while the two others are polarized 
on Poincare waves. This is done in Section [3] by proving some necessary conditions for the 
existence of those propagators. The general strategy is the following: 

(1) Consider the system A{x2,£D)U = irU. Take the Fourier transform in xi, which is 
possible since the equation is translation invariant in xi. Then extract from this system 
(by linear combinations and substitutions) a linear equation on one component n^ of U, 
of the type /i(x2,e-D2;^i,e,r)nfc = 0. 

(2) The symbolic equation corresponding to the PDE writes h{x2,S,2',Ci,£,'T) = 0. It has 
three roots (with respect to r) , t±{x2,C]£) (Poincare roots) and Tji{x2,C;£) (Rossby 
root). We find T±(a;2,^;e) = t±{x2,C) + 0{e) and tr{x2,C;£) = £Tr{x2,C) + O(e^). 

(3) Those roots are not necessarily symbols. To guarantee these are indeed symbols one 
needs a microlocalization. Given a compact set and a truncation on that compact 
set Xi one can construct three operators T^ (via a general theorem, stated and proved 
in an abstract way in Theorem[3]of Appendix [A|) whose principal symbols are precisely 
the Rossby and Poincare symbols t± and r/j. 

2.3.2. Trapping of Rossby waves. Section U] is devoted to the study of the Rossby propaga- 
tor To, and in particular to the proof of result ^ in Theorem[TJ For the time scale considered, 
it is easy to see that the energy propagates according to the trajectories of the semiclassical 
Rossby hamiltonian r^j. The first step of the analysis therefore consists in studying the dy- 
namical system giving rise to those trajectories. It turns out that the trajectories are always 
bounded in the X2 direction. One is therefore reduced to studying the trajectories in the xi 
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variable and in identifying the set A of initial data in the cotangent space giving rise to trap- 
ping in xi. One then checks that A is of codimension one under some additional assumptions 
on u, and the last step of the study consists in studying more precisely the trajectories in 
some specific situations, in particular in the case of the betaplane approximation. 



2.3.3. Dispersion of Poincare waves. The next step of our analysis of wave propagation con- 
sists in proving, in Section [5l that Poincare waves propagate so fast that they exit from any 
(bounded) domain of observation on the time scale that we consider, which proves result ([2]) 
in Theorem [U Note that, because of the very long time scaling, usual tools of semiclassical 
analysis cannot be applied for the Poincare waves : we actually need deeper arguments such 
as the Bohr-Sommerfeld quantization to conclude. 



2.3.4. A diagonalization result. Once the Rossby and Poincare propagators have been well un- 
derstood, we can retrace the steps followed in Section [3] to prove that the necessary conditions 
on the scalar propagators are sufficient. The difficulty is that the operators 11 enabling one to 
go from the original system to the scalar equations and back (computing an approximate left 
inverse Q of 11 at the order 0(e°°)) are only continuous on microlocalized functions; moreover 
the scalar propagators T^ are themselves only defined on microlocalized functions. So that 
requires understanding the persistence of the microlocalization of the solutions to the scalar 
equations. That is achieved in the two previous sections, where it is proved that if the initial 
data is conveniently microlocalized, then for any time t > one can find a compact set K 
and one can construct T^ as in Section [3l so that the solution to the scalar equations with 
propagators T^ is microlocalized in K (actually for Poincare modes the microlocalization is in 
the variables (a;2,Ci5C2) only, which is enough for our purpose). This enables us in Section [6] 
to conclude rather easily by computing explicitly the matrix principal symbols of 11 and Q. 



2.3.5. The analysis of the nonlinear equation. Section [7] is devoted to the proof that the 
solution of the nonlinear equation remains close to that of the linear equation. The method of 
proof consists first in proving the wellposedness of the nonlinear equation on a uniform time 
interval, by using semi-classical weighted Sobolev type spaces, whose additional feature is to be 
well adapted to the penalization operator A{x2,£D): one therefore constructs a matrix- valued 
pseudo-differential operator which approximately commutes with A{x2,£D). The convergence 
of C/e — V^ to zero relies on a standard L^ energy estimate on [7^ — 14, and Gronwall's lemma. 



2.3.6. Two appendixes. In AppendixElone can find the statement and the proof of a general 
theorem, used in Section [3l allowing to associate to a linear evolution PDE a number of opera- 
tors describing the dynamics of the equation; those operators are constructed by writing down 
the symbolic equation associated to the PDE and in quantizing the roots of that polynomial 
(in the time derivative) . Appendix [B] finally collects a number of prerequisites on microlocal 
and semiclassical analysis, that are used throughout the paper. 



semiclassical and spectral analysis of oceanic waves 
3. Reduction to scalar propagators 

3.1. Introduction. Let us first recall that the propagator 

(eu ■ eV edi £82 

edi eu • eV -h + e^u\ 
ed2 b eu- eV 

can, in the particular case when n = and 6(x2) = /3x2, be diagonalized without any error 
term (in particular for any finite e), using a Fourier basis (exp(-xi^i))^^gR, in xi and a Hermite 
basis ('0^(j;2))n6N in X2. Precisely, the following statement is proved in [13]. 

Proposition 3.1 (Gallagher & Saint-Raymond, [13]). For all {^i,n,j) G R x N x { — 1,0, 1}, 
denote by T(^i,n, j) the three roots (in increasing order in j) of 

(3.1) T^ - i^l + ^e{2n + 1))t + e^^i = 0. 

Then there exists a complete family of L^(R x R, R'^) of pseudo-eigenvectors (^5 ^ ■) of the 
operator A{x,eD) (where u = and 6(^2) = I3x2): 

(3.2) V(ei,n,i) G R X N X {-1,0,1}, A{x,eD)^^^^j = iT{^i,n,j)^^^^^^ 

where ^f ^ • can be computed in terms of the n-th Hermite function ^^(x2) and its derivatives. 

In other words, the three scalar propagators (numbered by j) can be obtained from the 
symbolic equation (j3.ip remarking that (3e{2n + 1) is the quantization of the harmonic oscil- 
lator — e^c^l + P'^^'q.- It is proved in [13] that as ^1 and n go to infinity 



r(a,n,±)~±We? + /?e(2n + l), and T(ei,n,0) 



e/36 



ii + /5e(2n + 1) 



We are interested here in deriving a symbolic equation similar to (j3.ip for a general zonal 
current u = (mi(x2),0) and Coriolis parameter b = 6(3:2). The difficulty comes from the fact 
that the propagation of waves is governed by a matrix of differential operators with non- 
constant coefficients, the diagonalization of which is not a standard computation. Of course, 
in the semiclassical limit e — t- 0, we expect to get a good approximation of the propagation at 
leading order by considering the matrix of principal symbols 

/O i^i 
(3.3) Aq{x,0-= Uii 

\«6 b{x2) 

and by computing the scalar propagators associated to each eigenvalue 




(3.4) r±(x2,6,6) :=±^ef + ^i + &'(^2) 

and 0. The eigenvalue corresponds to the Rossby modes, whereas the two 0(1) eigenval- 
ues ±Y^^^ + Ci + b'^{x2) are the Poincare modes. Nevertheless, this approximation is relevant 
only for times of order O(^), and we are interested here in much longer times, of order 0{\). 
This means that we need to compute the next order of the expansion of the eigenvalue 0. Once 
that is done, we need to quantify these symbol eigenvalues to deduce scalar propagators. 
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3.2. The Rossby modes. Finding precisely the Rossby modes (up to an 0{e^) error) requires 
more intricate calculations than merely diagonalizing the matrix of principal symbols Ao{x, ^) 
given in (|3.3p . So let ir be an eigenvalue of the propagator, assumed to be of the form r = 
etr + O(e^). Explicit computations lead to the subsystem 

Ci eui^i-Tj\uiJ \-i{b-e'^u[)u2^ 

and defining ar{x2,^i) := £Ui^i — r and Pr{x2,Ci) = —^i + ctr to the scalar equation 

(3.5) [ed2P~^^{iared2 + i^iib - e^u[)) - bp^{i^i£d2 + iar{b - e^u'^)) + ior) U2 = 

where .^i is the Fourier variable corresponding to xi/e. From now on we assume that ^i is 
fixed, and is bounded away from zero (recalling Assumption ()2.3p ). Note that equation (j3.5p 
makes sense because we are assuming here that r = etr + O(e^), so p^ is well defined. That 

would not be the case for the Poincare modes (where r = ±\/|^P + b'^{x2) + 0{e)) so we shall 
use another subsystem in the next paragraph to deal with the Poincare operators. 

In order to derive a symbolic equation associated with the differential equation (13. Sh . we shall 
proceed by transforming (|3.5p into a differential equation which is the left quantization (in the 
sense recalled in Appendix[B]) of a symbol, polynomial in e and r. This leads to a differential 
equation of the type 

(3.6) 'P2{x2;e,T){ed2)'^U2 +Vi{x2;e,T)ed2U2 + 'Po{x2;e,T)u2 = 0, 

where each Vj{x2',s,t) is a smooth function in X2, and has polynomial dependence in e and 
in r (precisely of degree at most 5 in e and r). This generalizes (|3.ip : one can compute in 
particular, using the fact that ed2ar = O(e^), that 

V2{x2;£,t) = ipRar + 0{e'^), 'Pi{x2;e,T) = 0{e'^), and 

Po(x2; e, t) = ipR {Cieb' - (6^ + ^f)ar) + 0{e^). 
The differential operator appearing on the left-hand side of (j3.6p is the left quantization of 
the following symbol: 

(3.7) h{x2,S.2;£,T) := -V2{x2;e,T)S.i + i'Pi{x2;e,T)S.2+'Po{x2;e,T) 
which belongs for each r to 5*2 (^r) for some function Qt of the type 

gAx2,^2) = {l + r'){l + b\x2)+e2) 

recalling that ui and u'l are bounded from above, as well as Assumption ()2.2p . 

Now we recall that h{x2,S,2',£, t) is a polynomial of degree 5 in r hence has five roots, among 
which 2 are actually spurious: these are of the form =b^i + 0(e), and they appear because we 
have multiplied the equation by a polynomial in r which cancels at the point it^i at first order 
in e. An easy computation allows to obtain the two other 0(1) roots, which are precisely the 
Poincare roots ±\/|^P + b'^{x2) + 0(e), and an asymptotic expansion allows also easily to 
derive the Rossby 0{e) root: one finds 

TR{x2,ii,i2]e) := erR(x2,Ci,6;e) + 0(e^), where 

(3.8) rij(x2,ei,6) := ^^^^-^^ + ^iMx2)- 

Now that the root tr has been computed, our next task is to prove the existence of the Rossby 
propagator Tr = eTr whose principal symbol is precisely er/j. Actually this result is a direct 
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consequence of Theorem [3] stated and proved in AppendixO with the notation of Theorem [3l 
one has u = I and drho{x,^,0) = —iCiiCi +?2 +b'^{x2))- In the foUowing statement, the time 
variable s is defined as s = t/e^. 

Proposition 3.2 (The Rossby propagator). Let tr he the symbol defined in Ii3.^) . Then 
for any compact set K, satisfying Assumptions i2.3^2.4\ ) there exists a formally self-adjoint 
pseudo-differential operator Tr of principal symbol tr such that if (pR is microlocalized in fC 
and solves 

(3.9) dsipR = ieTRipR, 

then 

'P2{x2;e,ds){ed2f^R + 'Pi{x2;e,ds){ed2)(pR + Vo{x2;e,ds)ipR = 0{e°°). 

Definition 3.3 (The Rossby operator). We shall call Hr the Rossby operator defined by 

/ P^^{i{-ie'^uidi - efR)ed2 + edi{b - e'^u'-^)) 
nR:= -p-^^-e^di2 + i{-ie^uidi-eTR){b-e^u'^)) 
\ Id 

where Pr := e'^df — {ie'^uidi + eTrY . 

Remark 3.4. (1) Notice that Hr is well defined since the principal symbol of Pr is 
bounded from below (see Appendix W\) . 

(2) Proposition \4-6\ shows that if ipRU-Q is microlocalized in a compact set /Cq satisfy- 
ing 1^2. 4\ ), then the solution to i3. 9\) is microlocalized for all t > in a compact set ICf 

(3) The above computations allow to formally recover the original shallow-water equation, 
up to 0(e°°). Indeed retracing the steps which enabled us above to derive equation /i3.6\) 
shows that if ip^ is a smooth function conveniently microlocalized, and if ip solves 

edtip = if Rip, 99|i=o = 9?°, 

then the vector field U := URip satisfies \2.1\) up to 0(e°°). This property will be made 
rigorous in Section\^ 

3.3. The Poincare modes. In this paragraph we shah fohow the method used above in 
the case of Rossby modes to infer Poincare propagators T± and operators II-i-. Actually 
one cannot use precisely the same method since the symbol (eui^i — r)^ — ^^ may vanish 
when T = T± + 0{£). So we shall instead consider the subsystem 

'euiii - T -ied2 \ [ p\ ^ f-S,iui 
-ied2 eui^i - tJ \U2j \ ibui 

and the scalar equation 

(3.10) (CiPpH-(^tCi - ed2h) + a^ + {b - e^u'i)p-p^{ed2 - a^6))ni = 



where as before ar{x2,ii) := £Ui^i — r, and ppix2,^i) := e^5| + a^. We notice here that a~^ 
is well defined when t = t± since t± is bounded away from zero by the assumption on ^i, 
and the same goes for r^ — ^| so pj, is also well defined. Then it remains to follow the steps 
of Paragraph 13.21 to obtain a new scalar PDE of the same type as (j3.6p . as well as a symbol 
equation of the type (j3.7p : 

(3.11) h{x2,C2;£,T) := -V2{x2;e,T)il + iVi{x2;e,T)i2+'PQ{x2;e,T). 
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Of course t± are roots of that equation up to 0(e), and etr is a root up to O(e^). Then the 
appHcation of Theorem [3] impUes a similar result to Proposition 13.2^ noticing that with the 
notation of Theorem^ v = and drho{x,^,T±) = —2{^f + ^| + b'^{x2))- 

Proposition 3.5 (The Poincare propagator). Let t± be the symbol defined in ( [3.^[ j. Consider 
a compact set ICp C R* x T*R. Then there exists formally self-adjoint pseudo-differential 
operators T± of principal symbols t± such that if (p± solves ds'^± = iT±ip±, and ip± is mi- 
crolocalized in R x /Cp, then 

'P2{x2;e, ds){ed2f^± + A(X2; e, ds){ed2)^± + Mx2; e, ds)ip± = 0(e°°). 

Definition 3.6 (The Poincare operator). We shall call Il± the Poincare operator defined by 

I P^^{i{-ie^uidi-T±)edi+£d2h) \ 
n± := Id , where P± := {ie^uidi +T±f + e'^dl 

\ P^^{i{-ie^uidi - T±)b - e^ieaa) / 

Remark 3.7. As in Remark \3.4l one sees formally that ifip^ is a smooth function conveniently 
microlocalized, and if (p solves 

e^dtip = iT±ip, ip\t=o = ^^, 

then the vector field U := Il±ip satisfies h2.1\) up to 0(e°°). This property will be made rigorous 
in Section\^ 

4. Study of the Rossby waves 

4.1. The dynamical system. For the time scale considered here, the propagation of energy 
by Rossby waves is given by the transport equation (see Appendix [B|l 

dtf + {rR,f] = Q 
where tr is the principal symbol of the Rossby mode computed in (jS.Sp : 

(4.1) y,(j„,„5,)._^J|2)|l_ + a,,,,)5,, 

As tr is a smooth function of (x2, Ci) C2)) the energy is propagated along the bicharacteristics, 
i.e. along the integral curves of the following system of ODEs: 

Since the condition (j2.3p avoids the set {^1 = 0}, we can suppose that ^j* ^ 0. Moreover 
since tr does not depend on xi, we find ^* = '^i- The ODE to be studied is therefore 



(4.2) 



^.^,^(^.)^6'(4)(-ef + ef + 6^(4)) 



■ t _ -2b-(x*,)6^2 
' iil + it + b^iA)? 



{il+e2+b^{A)? 



ui{x\)ii + 



2h{xi)b'{x\fii b"{xi)C, 



Due to the assumptions on ui and on b, the map (x,^) 1— ;• (V^r^, —'VxTr){x,C) is bounded, 
so the integral curves are globally defined in time. The strategy to study their qualitative 



SEMICLASSICAL AND SPECTRAL ANALYSIS OF OCEANIC WAVES 13 

behaviours is to first (in Section 1^2]) consider the motion in the reduced phase space (a;2, C2) £ 
R^ and then (in Section [4. 3p to study the motion in the xi direction. 

4.2. Trajectories in the reduced (a;2)'^2) phase space. In this section we study the tra- 
jectories in the reduced (2:2,^2) phase space. We shall denote ^1 := ^5*. 

4.2.1. Energy surfaces. Since the Hamiltonian tr and ^1 are conserved along any trajectory, 
trajectories are submanifolds of 

£t,^i ■= {(a;2,6) e R-^;T7?(6,a;2,6) =r}. 

In the following we shall note for any energy r and any ^1 G R* 

K,.(X2):= ''['f' -il-b\..), 



SO that if P := {x2/Vt,^i{x2) > O}, then Er,^^ = \{^2-,^JyT,(,i{x2)), X2 G V\. Note 
that Vr^iiix^) becomes singular if Xg reaches a point X2 such that r = ui{x2)ii. 
Proposition 4.1. The projection of £r,^i on the X2-axis is bounded. 



Proof. We recall that on Sr,^^ we have 



b'{x2)i 



^ +Ul{x2)il =T. 



e,+^i+b^{x2) 

Suppose the trajectory in X2 is not bounded, then in particular it escapes the support of ui. 
In the case when r 7^ 0, letting \x2\ go to infinity yields a contradiction due to the assumptions 
on b. In the case r = 0, for X2 out of the support of ui we have 

e,+e2+b^ix2) 

and the only possibility is for X2 to be fixed on a zero point of b' (hence in particular does not 
go to infinity). CH 

Proposition 14.11 shows that to prove that some trajectories are trapped in physical space, it 
suffices to study their behaviour in the xi direction. However before doing so, let us prepare 
that study by classifying the trajectories in the reduced phase space. Up to a change of 
parameter, namely expressing time i as a function of X2 

dt = ± ^'^^'^^\ dx2 , 

(r -^i'Ui(x2))V^^.5i(^2) 

which is justified locally, and will give the convenient global behaviour by suitable gluing, we 
are brought back to the study of the hamiltonian system ^1 ~ ^,§1(3^2) describing the motion 
of a particle in the potential — VV,gi. 

For smooth potentials V such that V{x) — )• —00 as \x\ — )• 00, the possible behaviours of 
such a system are well-known and the trajectories are usually classified as follows (see for 
instance [U [21 [HJ [20l [22]): periodic orbits, fixed points, homoclinic and heteroclinic orbits 
connecting unstable fixed points. Here the situation is more complex insofar as Vt,^i admits 
singularities. We shall classify the trajectories according to their motion in the X2 variable. 
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4.2.2. Periodic trajectories. These correspond to the case when there exists [xmin,Xmax] in R 
with Xmin 7^ Xmax , Containing X2 such that 

• yr,5i has no singularity and does not vanish on ]xmin,Xmax[', 

• the points Xmin and x^ax are reached in finite time. 

The extremal points Xmin and Xmax are then turning points, meaning that the motion is 
periodic. The fact that Xmin and Xmax are reached in finite time is equivalent to V^c_^ (xmin) 7^ 
and V^^^_^{xmax) + 0. Indeed if V'^^^^{xmax) = (resp. V^^^_^{xmin) = 0), then {xmax,0) (resp 
{xmin,0)) is a fixed point, which contradicts the uniqueness given by the Cauchy-Lipschitz 
theorem. And conversely, if V^ f_^{xmax) 7^ (resp. V^ f_^{xmin) 7^ 0), an asymptotic expansion 
in the vicinity of x^ax (resp. Xmm) shows that the extremal point is reached in finite time. 

Definition 4.2. We will denote by V the subset of the phase space T*R^ consisting of initial 
data corresponding to periodic motions along X2- 

A rather simple continuity argument allows to prove that V is an open subset of R x R* x R^. 
Denote indeed by (ai?,'^!, 3;2,'^2) ^^y point of P, and by Xmin and Xmax the extremal points 
of the corresponding (periodic) trajectory along X2. As VI' - (xmax) 7^ and Vt^^^ is a smooth 
function of r and ^1 outside from the closed subset of singularity points, the implicit function 
theorem gives the existence of a neighborhood of (^1 , f ) such that there exists a unique x^ax 
which satisfies Vr^^-^ixmax) = 0. Furthermore Xmax depends continuously on r and ^1, in 
particular V^c {xmax) 7^ 0. Using the same arguments to build a suitable Xmim we finally 
obtain that there exists a neighborhood of {x^,^i,X2,$.2) ^^ which the motion along X2 is a 
non degenerate periodic motion. Moreover, Xmin, Xmax and also S,min, Cmax and the period T 
depend continuously on the initial data. 

Fixed points correspond to the degenerate case when Xmin = Xmax = 2^2 > which implies 
that either ^2 = or b'{x2) = 0. The latter case is completely characterized by the condi- 
tion b'{x2) = 0, so let us focus on the case when b'{x2) 7^ 0. Fixed points correspond then to 
local extrema of Vr,5i • They can be either stable or unstable depending on the sign of V^'c ■ 

Stable fixed points are obtained as a limit of periodic orbits when the period T — )■ 0, whereas 
unstable fixed points are obtained in the limit T — )• 00 as explained below. 

Stopping trajectories belong to the same energy surfaces as unstable fixed points and reach 
some unstable fixed point in infinite time : they correspond to the case when there exists an 
interval [xmin,Xmax] of R containing X2 such that 

• Vr^^-i has no singularity and does not vanish on ]xmin,Xmax[] 

• Xmin and Xmax are either zeros or singularities of Vr^^-^ ; 

• ast^oo,xi^xf such that V^^^^{xf) = or b'{xf) = 0. 

They are in some sense also a degenerate version of periodic trajectories since for arbitrarily 
close initial data, one can obtain periodic orbits. 

Definition 4.3. We will denote by 6V the subset of the phase space T*R^ consisting in initial 
data corresponding to fixed points and stopping motions along X2 ■ 
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Using the characterization of the energy surfaces which carry such pathological motions, we 
can prove that 6V is a codimension 1 subset of the phase space. We first consider the energy 
surfaces containing a fixed point X2 such that &'(x2) = 0. We then have r = £,iui[x2)- The 
corresponding set of initial data 

{{xl,ii,xl,il) /TR{ii,xl,il) =6^1(^2)} is of codimension 1. 

As we have assumed that h has only a finite number of critical points, the union of these 
sets is still of codimension 1. We then consider the energy surfaces containing a fixed point 
(xi,^i,X2;0) such that V- ? (X2) = V- x (^2) = 0. We then have 

-tf.fe)+ ^,''fe"''fe' -.""fe' =0. 

(«? + »(i2))-' (f + (,= (%) 

For each X2, there are at most two values of ^f such that the previous quantity vanishes. We 
therefore deduce that 



„o c ^0 tO^ /~^rt ^0 cO\ _ b'{x2)^i 

e? + fe2(x2) 



(x?,gi,x°,g°)/Tfi(6,x°,e°)= ,2 'A N +^1(^2)6 and 



_,^, ^ , 2bix2)b'{x2) b"{i2) n\- f H- • 1 

— Ui UC2 H — ^ ^ = U > IS 01 codimension 1. 

iei+bHx2)r ei+bHx2) -< 

It consists indeed of at most eight manifolds, each one of them parametrized by the real 
parameter X2- In the sequel, we shall avoid these pathological motions assuming that the 
initial data is microlocalized outside dV. 

4.2.3. Asymptotic trajectories. These correspond to the case when there exists an inter- 
val [xmin,Xmax] of R Containing X2 such that 

• Vr^^-^ has no singularity and does not vanish on ]xmin,Xmax[', 

• Xmin and Xmax are either zeros or singularities of Vr^(^ ; 

• as t — )• 00, X2 — )• ^2° where ^2° G {xmin,Xmax} is a pole of multiplicity 1 of V^.^^. For 
the sake of simplicity, we further impose that 6'(x2°) 7^ 0. 

This situation therefore corresponds to a motion which is not periodic. 

Depending on the sign of ^1, an asymptotic trajectory will either encounter a turning point 
and then converge asymptotically to the singular point, or converge monotonically to the 
limiting point. As X2° is such that ui(x2°) = t/^i, one has 

• either r 7^ and x^ belongs to the support of ui, 

• or T = and Ui{x2) = —b'{x2)/{^f + ^2 + ^^(^2)) / 0' meaning that Xg belongs to 
the support of ui. Therefore, X2° is either min{y > X2 / ui{y) = 0} or max{y < 
X2 / ui{y) = 0}, in particular x^ belongs to the support of ui. 

Definition 4.4. We will denote by A the subset of the phase space consisting of initial data 
corresponding to asymptotic motions along X2- 

The same kind of arguments as in the previous paragraph allow to prove that A is an open 
subset of the phase space R x R* x R^. Consider indeed some (x5,Cij^2'?2) ^ -^^ and the 
corresponding asymptotic point x^. As x^ is a pole of multiplicity 1 of V^ ; , the implicit 
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function theorem shows that V^^ admits locally a unique zero, which depends continuously 
on T and ^i. Using further the continuity of the possible turning point, we get that A contains 
a neighborhood of (X]',^i,i2'^2)- Moreover, we can obtain bounds on the expansion (with 
respect to time) of any compact subset of A. Here we will focus on the growth of ^2; s-^d 
proves that it depends continuously on the initial data in A. Without loss of generality, we 
can consider the case when the asymptotic point is Xmax- Then we recall that 

lim ^2 = oo and lim X2 = x^, with ^1^1(^2°) = t- 
As X tends to ^2°, we have (recalling that b'{x2^) / 0) 

This implies that 

\tt\2 "1^2 ) / oo\-l T -t ^1,1/ oo\c ^2 0^1 '''"2 ) Si / oo \3/2 

S2 ~ , , ^x (x — Xo ) and Xo ~ — 26 (xo jciT-nr ~ 2 — — —7- — (xo — x) ' . 

By integration, we get 

X*2~xf=+C7ii-2, e2~C'2t 

where Ci and C2 depend continuously on x^, and consequently on the initial data. 

Singular trajectories are a degenerate version of the asymptotic trajectories above: they cor- 
respond to the case when there exists an interval {xmin-,Xmax\ of R containing X2 such that 

• Vt,s,i ^^^ '^o singularity and does not vanish on ]xmin,Xmax[', 

• Xmin and Xmax are either zeros or singularities of Vr', 

• as t — )• 00, X2 — )• x^ where x^ E {xmin-,Xmax} is either a singularity of order greater 
than 1 of VT,i-i_ , or a singularity which is also a zero of b' . 

As previously, it is easy to check that x^ is necessarily in the support of ui. Furthermore, 
we have it'^(x^) = 0. Indeed, in the vicinity of x^, we have 

,. , , b'ixf)Ci 



a{xf) - u{y) 

which has a singularity of order greater than 1, or a singularity which is also a zero of h' if 
and only if u'(x^) = 0. 

Definition 4.5. We will denote by 5 A the subset of the phase space consisting of initial data 
corresponding to singular motions along X2- 

The previous condition tii(x2°) = shows that 6 A is included in the union of energy surfaces 

£t,£_i with r = uiiij) for some y G supp{ui) such that n'^(y) = . 

As u'l has only a finite number of zeros in the support of ui, this implies that 6 A is a 
CO dimension 1 subset of the phase space. 

Gathering all the previous results together, we obtain the following 
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Proposition 4.6. The phase space R x R* x R^ admits the following partition 

RxR* xK^ = VUAUT, 

where V and A are the open sets of initial data giving rise respectively to periodic motions 
and asymptotic motions along X2, and S := 6V U 6 A is the codimension 1 set of initial data 
giving rise to pathological motions along X2. 

For any compact set K C AUV, and for any time T > 0, we further have a uniform bound 
on the image of K by the flow up to time T. 

Proof. The first statement just tells us that all trajectories belong to one of the four categories 
described above. Indeed, for any initial data {xi,^i,X2,C2)j o^^^ has 

^T,^i{x2) = (C2) ^ 0) s-^d Vt,^i{x2) — ^ —CO as \x2\ — )■ 00, 

so that there exists an interval [xmin,Xmax] of R containing X2 such that V^^^^ has no singularity 
and does not vanish on ]xmin,Xmax[, and Xmin and Xmax are either zeros or singularities 
of V^,^i • The second statement is then a simple corollary of the continuity results established 
on V and A. D 

4.3. Analysis of the trajectories in the xi direction: trapping phenomenon. Propo- 
sition [54] states that the trajectories are always bounded in the X2 variable, so it remains to 
study the xi variable. Our aim is to find a set A C R x R* x R^ such that any initial 
data {xi,^i,X2,C2) in A gives rise to a trapped trajectory, meaning that 

xfds is uniformly bounded for t G R^ . 


4.3.1. The criterion of capture. Let us prove the following result. 

Proposition 4.7. A necessary and sufficient condition for a trajectory with initial data in AU 

V to be trapped is 

1 /■* 
lim - / xfds = 0, 

where T denotes the (finite) period of the motion along X2 in the periodic case, and T = +00 
in the asymptotic case. 

Proof. We will study separately the different situations described in the previous section, 
namely the case of periodic and asymptotic trajectories in (3:2, ^2)- 

• In the case of a periodic motion in (x2,^2) of period T > 0, the function x\ is also periodic, 
with the same period. Writing 

x\ = ^? + / ( *i " ^ / ^1' '^A "^"^ + ^ / ""1 "^^ 

we see that depending on the average of x\ over [0,T], x\ is either a periodic function, or 
the sum of a periodic function and a linear function. It follows that trapped trajectories 
are characterized by the criterion /„ x\dt = 0. Note that, depending on the period T, the 
trajectory can explore a domain in xi the size of which may be very large. Nevertheless the 
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continuity statement in Proposition 14.61 shows that we have a uniform bound on this size on 
any compact subset of V. 

• For asymptotic motions, we need to check that 

x\ — ui{x^)t is integrable at infinity. 
We have indeed 

x] = uiixi) + ^ ^^^ ^ „ ^ \_j^ 

' " {il + it + h^iA)? 

which, together with the asymptotic expansions of x^ and ^2 obtained in the previous section, 
impHes that 

x\ = ui{x^) + 0{r^). 

It is then clear that the trajectory is trapped if and only if 

ui{xf)= lim - / ±lds = {). 

t^cc, t Jq 

Using again the continuity statement in Proposition 14.61 we also get a uniform bound on the 
size of the time evolution of any compact subset of A. D 

Remark 4.8. Note that, in the case of a singular asymptotic motion along X2, the criterion 
of capture is equivalent to u{x^) = t = 0. In the case of periodic motions along X2, this 
criterion - even more complicated - can also be expressed in terms of the initial data. The 
map x\ is indeed a smooth bijection from a time interval ]tQ,tQ -\- ^[ to ]xmin,Xmax[, so 

T f-to+T/2 rXmax 

±\ dt = 2 x{dt = 2 xi\.^^ (i2)ri=.. dx2 . 



The trajectory is therefore trapped if and only if 

(4.3) /-[j^_^«;)]^„, ,,)-./..,.„ 

Jx,r,i„ L 2((r/^i) -ui{y)) J 

(notice that since we are only looking for a criterion for the function to vanish, we can re- 
place ^2 by Y^Vr,5i without discussing the sign). Such a formula would be useful to investigate 
numerically the initial data giving rise to trapped trajectories. We shall also use it when 
investigating in more detail the case of the betaplane approximation, in Paragraph I 



4.3.2. Exhibiting a subset of A of codimension 1. Let us prove the following proposition. 

Proposition 4.9. Let ui G (7^(R) be a function which is not identically positive, with some 
zero of finite multiplicity. Then the set Aging consisting of data in A giving rise to singular 
and trapped trajectories, is nonempty. It contains a submanifold 0/ R x R* x R^ which is of 
codimension 1. 

Proof. Without loss of generality, we can assume that there exist yi < y2 where ui vanishes, 
with u'i{y2) > and ui{y) < on ]yi,y2[- As we want to study trapped asymptotic trajecto- 
ries, we will restrict our attention to the case when r = 0, which is a necessary condition for 
asymptotic trajectories to be trapped. Extremal points of the trajectories are then defined in 
terms of the function 



^o,a(2/)=-P^-%)^-d 
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More precisely, if we introduce the auxiliary function 

Q{y) ■■= - ^-r^ - b{yf ^ yG]yi,y2[, 
ui{y) 

we obtain turning (or stopping) points y^ if piVs) = Ci: and singular points ys if lim p{y) = 
+00 or equivalently ui{ys) = 0. By definition of yi and 1/2, one has 

lim Q{y) = +00 and lim Q[y) = +00. 

Let us define N := max ( ; inf Q{y) ) G R+ . For ^1 such that ^f > N , we then define 

J/G]yi,J/2[ 

h{ii) := sup {y g] -oo,y2[; £'(y) < Ci } G]2/i,y2[- 

We therefore have that 

Vy G]/i(Ci))y2[) y is neither a turning point nor a singular point. 

As /i is a decreasing function on ] — 00, — \/]V], all ^1 g] — 00, — \/iV] except a countable number 
are continuity points. Choose then some .^i to be a continuity point of h and X2 £]^(Ci)) ?/2[- 
By continuity of /i, there exists a neighborhood V of (521^1) such that 

The set {(x?,^,^^, {q{xI) - ^i)^ ) ; (x?,^,^:") G R x T/} is a submanifold of R x R* x R^ 
having codimension 1. Furthermore, for any initial data in this set, we have i;2|t=o > and 
a simple connexity argument shows that x^ is an increasing function of time. In particular 
a^2 ~^ 2/2 as t — )■ 00. This proves Proposition 14.91 D 

4.4. Some examples in the betaplane approximation. We are concerned here with the 
betaplane approximation, that is when b{x2) = 13x2- One has therefore 

M^l,X2,^2) = Ui{^2)^l + ^2 ^ ^2 I p2^2 

and with the notation of Paragraph 14.21 



K,a(^2) = ^-e?-/3V 



r - 2- 



• In the absence of convection, one can characterize exactly the set A of initial data giving 
rise to trapped trajectories. One can notice that for r such that r^^ G ]0,/3,^j~ ], the energy 
surface 1?,-^^^ is simply the ellipse 

{(^2,e2)/ei + /3'^i = /3eir-^-ei}. 

Let us go through the previous analysis and study the trajectories in this situation. One 
notices that fixed points correspond to X2{t) = C2(i) = 0, with xi{t) = x^. There are no 
asymptotic trajectories (singular trajectories would correspond to r = 0, which is not possible 
here since ^1 7^ 0). Finally let us consider periodic trajectories. In order to get trapping one 
must check that if the energy level r and the frequency ^1 are fixed, we have 

(3(^-2^1) 

dx2 = 0. 



i 



^^-ef-/32xi 
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Let US choose for instance any (xi,.^i) G R x R* and define the energy level r = /3/(2^i). 
Then the integral is identically equal to zero hence the corresponding trajectory is trapped. 
It corresponds to x\ = x^, and to X2,i^2 satisfying for all times ^2 + f^'^^2 ~ ^i- 

Such an explicit characterization cannot be obtained if there is some convection, in particu- 
lar in physically relevant situations. Nevertheless, we are able to prove that, under suitable 
assumptions, A is not empty, and more precisely that it contains both singular trapped tra- 
jectories and periodic trapped trajectories. 

• Proposition l4.9l shows that under the assumption that ui G C^(R) is a function which is not 
identically positive, with some zero of finite multiplicity, A contains a subset of codimension 1 
of initial data giving rise to singular trapped trajectories. 

• Finally let us construct periodic trajectories. We suppose to simplify the notation that /3 = 1, 
and that ui has a local maximum at zero, with, say 

< ni(0) < 2/3 and <(0) < -6. 

Then we shall prove that the set Aper is nonempty and contains a submanifold of R x R* x R^ 
which is of codimension 1. Let r be an energy level. We notice that for r] small enough 
and y g] — 7], ry[ 

ui{y) <ni(0) < ^ 

so there is no singular point for 1/^,^1 in the interval ] — Ty,ry[. Now let us show that there are 
two turning points inside ] — rj,ri[ for r/ small enough; this will imply that there is a periodic 
trajectory of energy r. Define the function 

Then H'^ (0) = and an easy computation shows that if rj is chosen small enough, then for 

all (y,6) G] -^,^[x[1,+oo[, 

^«i(y) = Sf^-<'(y)>V2>0. 

^ i^t + y) 

Now let us choose r] such that sup Ui{y) < —3, and let us define the number 6 := 1 . 

]-ri,v[ ^ 

c -1 e 

We choose from now on ^1 > and r = t(^i) := ui(0)^i H so that H^^{0) = ^r— < 0. 

The function H^-^ is decreasing on ] — rj,0[ and increasing on the interval ]0,r/[; moreover 

H^,{±^) > if^^(O) + ^ = -(1 - 5) + ^ > ^ > 0. 

It follows that, for all ^1 G [1, -|-oo[, there are two points XminiCi) (in ] — ??,0[) and XmaxiCi) 
(in ]0,?7[) satisfying the requirements of periodic trajectories, in the sense of Section r4.2.2[ 

According to the criterion (j4.3p . we define the following function on [1, -|-oo[: 
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and let us prove it vanishes. We notice that recalhng that ui(0) < 2/3, 

r^maxil) p , \ 1 1 

G(l) > / [l-[M0) + 6)^\Vr,^,{yr'/'dy 



6<52 



since r] g]0, 1[. On the other hand, when ^i goes to +oo, we find that Xmm(Ci) ~^ and 
Xmaxi^i) —5- so that recalhng that t(^i) = ^1(0)^1 — 5/^i, 

£4 pXmax{£,l) 



G{Ci)--nm^ Vr,^,{y)-'/' dy, 



so hm G(^i) < 0. By construction, the function G is smooth so there is some .^i belong- 

ing to ]l,+oo[ such that G(^i) = 0. We have generically (in r) G'{^i) 7^ 0. This implies 
that Aper 7^ 0) where Aper is the subset of A giving rise to periodic trajectories, and that Aper 
contains a submanifold of R x R* x R^ which is of codimension 1 . 

5. Study of the Poincare waves 

5.1. The strategy. In this section we want to prove the dispersion property of the Poincare 
polarisation, namely result (2) of Theorem [TJ We recall that the principal symbols t± are 
of order one, so one needs to study diffractive-type propagation, on a time scale of the or- 
der l/e^. Computing the bicharacterstics, as in the previous section, is therefore not enough 
to understand the classical flow. We shall instead rely on a spectral argument to prove that 
Poincare waves do disperse, and escape from any compact set in the physical space. The first 
step consists in taking the Fourier transform in xi (recalling that the problem is invariant by 
translations in the xi direction). We also recall that the data is microlocalized on a compact 
set C such that ^1 is bounded away from zero. Since t± = =t\/^2 + ^^(^2) + Ci > functional 
calculus implies that one can find classical pseudo-differential operators -ff±(^i) of principal 
symbols ^| + 6^(x2) such that T± = ±-\/-ff±(?i) + Ci- Let us now call A^(^i) and (p±i^i; X2) 
the eigenvalues and eigenfunctions of H±{S,i)- The following proposition will be proved in the 
next paragraph. 

Proposition 5.1. Let (j) he an eigenfunction of H±{^i), microlocalized on an energy surface 
which interstects C. Then (j) and its associate eigenvalue A are C°° functions 0/^1. More- 
over -d^jX is bounded on compact sets in ^1. 

5.2. Proof of Theorem [1](2) assuming Proposition 15.11 Let us now carry out this pro- 
gram. We consider an initial data denoted ip^, microlocalized in C. One can take the Fourier 
transform in xi which gives 

/(x) = -^ U\i^,X2)e-^'^ dCi. 

Now let us consider a coherent state (in Fourier variables) at [q^p) (see Appendix[B|), that is: 

(5.1) ^qpiCi) ■■= ——re* - e 2. . 

(7re)4 
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After decomposition onto coherent states we get 

Or \ 1 1 /" -0/ N jiiii-LLl _ «i-p)' , , ,^ 

if [X) = / tp^{q,p,X2)e E e 2^ dqdpdt,i 

where (p^{q,p,X2) := (v3gp|(^'^(-,X2))^2 • We notice that the integral over p and q is, mod- 
ulo 0{e°°), on a compact domain due to the microlocalization assumption on (p^. Finally 
decomposing onto the eigenfunctions (^j.(^i;x2) gives 

'^ (-^) = ^S^: — :t >^ / '/' (9,P;«;,6,a;2)v'±(Ci;a;2)e ^ e 2. dqdpd^i 
V27re(vre)4 ^ J 

where ^ {q,p;k, S,i,X2) := ( (/7-|-(^i; •)|(^ ('ZjP)")) • Note that the dependence of ^° on ^1 is 



L2 

only through the eigenfunction ip^, so ^'^ depends smoothly on ^1, as stated in Proposition l5.ll 

The sum over k contains 0{e~^) terms, due to the fact that A^(^i) remains in a finite interval 
(this will be made more precise in the next section, see Remark 15 .Sp . Now it remains to 
propagate at time t/e'^ this initial data, which gives rise to the following expression: 



2TTe TTe)i .J 



dqdpd^ 



V2vre(^e)-. , 

The stationary phase lemma then gives that this integral is 0(e°°) except if there exists a 
stationary point, given by the conditions: 

ii = p and £{xi — (?) ± 1 = 0. 

The second condition gives 

1 



2y'Ai(ei) +ef(xi - 9) = T^ (2p + %Ai(ei)) t. 

Therefore, since p 7^ and the \\. 's are bounded, with d^^\\{^i) = 0(e), there is no critical 
point for x\ in a compact set. Proposition 15.11 therefore allows to apply the stationary phase 
lemma and to conclude the proof of result (2) of Theorem [TJ Notice that the (fixed) losses 
in £ (namely the negative powers of e appearing in the integrals and the number of fc's in the 
sum) are compensated by the fact that the result is 0(e°°); it is important at this point that 
as noticed above, the function Tp^ depends smoothly on ,^1. 

5.3. Proof of Proposition 15.11 The first step of the proof consists in using the theory 
of normal forms in order to reduce the problem to the study of functions of the harmonic 
oscillator (as in Bohr-Sommerfeld quantization, paying special attention to the dependence 
on ,^1). The second step then consists in checking that the eigenvectors and eigenvalues have 
the required dependence on ^1. In the following we shall only deal with Hj^ to simplify, and 
we shall write H := H^. The first step relies on the following lemma. 

Lemma 5.2. Suppose that S,i lies in a compact set away from zero. There is an elliptic 
Fourier Integral Operator U : 'C(R) — )■ 'D(R), independent of ^i, and a pseudodifferential 
operator V{^i) with C°° symbol in ^i such that microlocally in any compact set K C r*R one 
has {VUy = [VU)-^ and 

VUHiVUy' = f[-e^di + xhCi;e) + 0{en 
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where 

while /o is a smooth bijection on R+ and the fj 's are C°° functions of I and ^i . 

Proof. The proof consists in using techniques linked to the isochore Morse lemma (see [B], 
[30]). Actually we introduce a canonical change of variables (the corresponding operator 
being the FIO U) allowing to pass from the variables (3^25^2) to action-angle variables (see [1] 
for instance). Let us make this first step more precise : we recall that the action variables are 
given by 



I := (b C2dx2, 

where the integral is taken on a constant energy curve H = h. The angle variables are then 
given by solving 9 = diS{I,^2), where K is the hamiltonian in the new variables, which only 
depends on /, and where S is defined by 

dSj=constant = i2dX2- 

Note that this is a global change of variables. It is now well known (see [7] for instance) that 
such a canonical change of variables is associated with an FIO U (independent of ^1 since the 
principal symbol of H does not depend on ^1) such that 

(5.2) UHU'^ = f^i-e^dl + xl) + eFi 

where f^ is a smooth, global bijection on R+. In this case we can actually write (see [7]) the 
following formula for U: for any L? function ip 

1 f ■ S(x2,i2)-Vi2 

Uip{x2) = 3- / e* ^ ao{x2,C2,e)u{y)dydC2 

(27re)2 J 

where oq is constructed so that U is unitary (up to 0(e)). Using those new coordinates and 
the formula giving U as well as the formula for the principal symbol of the adjoint (which 
is here the inverse) given in Appendix [HI it is not difficult to show that (15. 2p holds. Once 
the function f^ is obtained we proceed by induction: we first look for a symbol qi such 
that Qi := Op^{qi) satisfies 

e^Qif/^fZ-ig-^Qi = /0.i(-e252 ^ ^2) ^ ^2^^^ 

In order to compute qi we notice that 



P*yi 



UHU^^er''^^ = UHU'^ + / dtHt dt 

Jo 



where Ht := e''^'^^UHU^^e~''^^\ One sees easily that 
dtHt = i[Qi,Ht] 

= e^*Qi [Qi, f\-e^dl + xl) + eFr]e-''^^ 

so the principal symbol of dtHt is therefore the Poisson bracket e{qi, / (^2 + ■^2)}- O'^^ then 
remarks that the equation 

{q,f{xl + e2)}=9 
has a solution if and only if g has zero mean value (in action-angle variables) : one has indeed 
necessarily 

dif 
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Note that if / and g are smooth, then so is q since dif > everywhere. It follows that one 
can find qi simply by solving 

{qufixl + e2)} = fi-7i 
where /^ is the average of /i 

7i{I) = ^ I ^ fi{Re{x2,C2);Ci) de 

where Rg denotes the rotation of angle 9. This implies that with this choice of qi and 
writing Qi := Op^{qi) one has 

One proceeds similary at all orders. D 

The next step consists in using that lemma to check that the eigenvalues and eigenfunc- 
tions enjoy the expected smoothness properties. Actually this is rather straightforward since 
if {ip, A) satisfy 

H^P = Xi; + 0{e°°) 

then ip := VUtp satisfies 

hence 

This implies that 



VUHiVUy^ij = AV' + 0{e°°) 



A = A„, = / f e(n + -);^i;ej and V' = V'n = — /^n ( — f= j 

where h^ is the n-th Hermite function. The conclusion follows recalling that /(/;^i;e) = 
/o(I) + ^ £■'/,•(/; ^i;e) (where /o does not depend on ^i) and that each eigenfunction is 
microlocalized in a compact set. 

Remark 5.3. Since we have a complete spectral description of T±, with discrete spectrum 
for each given ^i, it is obvious that if the initial data is microlocalized in ('^i,X2,^2); then the 
solution to the equation e'^dtif = iT±ip remains microlocalized in the set of energy surfaces 
containing (^^i, 2^21 '^2)- Notice also that there are 0(l/e) eigenvalues in a compact energy 
surface. 

6. DiAGONALIZATION 

In this section we shall prove that the scalar propagators defined in Section [3] correspond 
indeed to a diagonalization of the original linear system (j2.ip . Let us prove the following 
proposition. 

Proposition 6.1. Consider a compact set /C C R* x r*R. With the notation of Defini- 
tions {3^ and \3.(A the operator 11 := (n_ 11/^11+) maps continuously H^ functions, microlo- 
calized in R X /C and satisfying i2.3\) . onto H^ . Moreover it is left-invertible modulo e°° , and 
its left inverse Q (modulo e°° ) maps continuously those functions onto H^. 
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Remark 6.2. This proposition, along with Proposition \4.6\ and Reniark l 5. 3\ showing the prop- 
agation of the microlocal support (in (^i,X2,'^2)j of Rossby and Poincare modes, prove the first 
part of TheoremUl 

Proof. The main step consists in showing that 11 does have a left inverse, and in computing 
its principal symbol. The construction of a left inverse can be done symbolically as follows. 
We first compute the matrix-principal symbol Vq of 11. One gets 

/ -Ci VeTP-i^2b _ih_ eiV^TP-i^^ \ 

il + h^ ii il + b^ 

1 -^ 1 



Vo 



\ ei+b^ ei+b^ / 



This shows that 11 maps microlocalized functions in R x /C onto H^ . A simple computation 
shows that 

therefore Qo = 'Pq exists. Let us call Qo the matrix obtained by Weyl quantization (term 
by term) of Qo- By symbolic calculus we have that: 

QoU = Id + eh 

Let us call Qi = —IiQq , where Xi is the (matrix) principal symbol of Ii. Again by symbolic 
calculus we have that 

{Qo + eQi)n = ld + e^l2 
where Qi has principal symbol Qi. Defining now Q2 = —T2Q.Q we get 

(Qo + eQi + e^Q2) n = Id + e^/a 

where Q2 has principal symbol Q2, and so on. This allows to invert the matrix II up to 0(e°°), 
and the principal symbol of the (approximate) inverse matrix Q is given by Qq: we have 



^ / ibi2 - iWe+b^ il + b^ -ibve + i^ + 66 

^ ' \ibC2 + CiVe+^ 6' + 62 ib,/e+^ + CiC: 



2 



Note that each term of the expansion of the symbol of Q is a polynomial of increasing order 
in 6 but that is not a problem due to the microlocalization assumption. The operator Q 
therefore clearly maps continuously H^ onto itself. The proposition is proved. D 



7. Control of the nonlinear terms for a weak coupling 

We are now interested in describing the behaviour of our initial nonlinear system, which 
includes the effect of the convection by the unknown (ni,n2). We shall prove Theorem [2] in 
this section. We recall that the system reads as follows: 

(7.1) e'^dtU + A{x,eD^)U + e^+'^Si{U)ediU + e^+'^S2{U)ed2U = 0, r? > 
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with 



(7.2) 



Si{U) 



Ui 


2^0 








Ui 








S2iU) 



U2 





\uq 





U2 





^Uq 





U2 



and U = {uo,ui,U2)- The usual theory of symmetric hyperbohc systems provides the local 
existence of a solution to (|7.ip in H^(R?) for s > 2, on a time interval depending on e a priori. 
Because of the semiclassical framework (fixed by the form of the initial data), it is actually 
natural to rather consider e-derivatives. Moreover, as the derivative with respect to X2 does 
not commute with the singular perturbation, we expect even the semiclassical Sobolev norms 
to grow like exp (— ), and therefore the life span of the solutions to (|7.ip to be non uniform 
with respect to e. That is the reason why the W^ spaces were introduced in p.7|) . 



7.1. Propagation of regularity for the linear singular perturbation problem. Let 

us first remark that derivatives with respect to xi do commute with A{x,£Dx), so we can 
propagate as much regularity in xi as needed. Extending a result by Dutrifoy, Majda and 
Schochet [TU] obtained in the particular case when 6(2:2) = 13x2, we will actually prove that 
there is an operator of principal symbol (^| + 6^)Id which "almost commutes" with A{x, eDx) 
in the semiclassical regime. 

• The first step, as in [10], is to perform the following orthogonal change of variable 

'uo + ui no - Ul 



U :-- 



-,U2 



V2 ' V2 
in order to produce the generalized creation and annihilation operators 

L±:=-^{ed2Tb). 

The system ()7.ip can indeed be rewritten 

e'^dtU + A{x,eax)U + e^Si{U)ediU + e^S2{U)ed2U = 



with 



L^ + ^n'A 

A{x,eDx):=\ euiedi-edi L_ - ^u[ 




L, 



and 



Si{U) :-- 



/ 3tto - -Ul 

2V2 



V 











uq — 3ui 







Up -Ul 

V2 / 



S2{U) :-- 



V2 
euiedi ) 



U2 







U2 



\ 



Uq + Ul Uq + Ul 



Up + Ui \ 

4 
Uq + Ul 

4 

U2 



Next, remarking that [e 82 — b , £82 ± 6] = ±2eb'{ed2 ± 6) it e b" , we introduce the operator 

^e^di - 62 + 2e6' 



D, 








£2Q2 _ 52 _ 2eb' 








e^dl - 62 
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We notice that D^ is a scalar operator at leading order. Moreover one can compute the 
commutator [Ds,A{x,eDx)]: we find 



[De,A{x,eD^)] 



^'£^di,eui]£di ^{[e'^di,u[] + 2eb'u[ - h") ^ 

v2 






^ ^b" -j=h" [e^dleu,]ed, j 



• At this stage, we have proved that 

(7.3) \D,,A[x,eD,)\ = 0{e^{ld - e^df - D^)) 

meaning that the commutator [D£,A{x,eDx)] is of order O(e^) with respect to the elliptic 
operator Id — e^df — D^. That implies that the regularity of the solution to the linear equation 

s'^dtV + A{x,eDx)V = 

can be controlled by an application of Gronwall's lemma: one has 

e^\\{Id-e^dj-D,)Vit)\\l, < emid-e^df-D,)Vo\\l2+Ce^ f \\{Id-e^dj-D,)V{s)\\l,ds , 

Jo 

where C depends on the VF^'°° norms of ui and b, so 

Wild - e^dl - D,)V(t)f^, < CWild - e^dl - B,)V^f^, e^*. 

7.2. Uniform a priori estimates for the nonlinear equation. Since the extended har- 
monic oscillator controls two derivatives in x^^ we get a control on the Lipschitz norm of \J 
of the type 

(7.4) \\ed,V\\i^^ < ^{\\D^,U\\l2 + \\e^dfUh2 + \\Uh2) . 

As Dg is a scalar differential operator at leading order in e, the antisymmetry of the higher 
order nonlinear term is preserved. More precisely, we have, using the Leibniz formula, 

e^dtDeil + A{x, eDx)DeU + e^Si{U)ediDeU + e^S2{U)ed2DeU 

= -[De,A{x,eDx)]U - e''S2{U)[D,,ed2]U - e^[D,,Sj{U)]edjU 
as well as 

e^dtDltj + A{x, eD.x)DlU + e^Si{U)ediD^^U + e^S2{U)ed2D'^U = -[D^,A{x, eD^)]U 
-e^S2iU)[De,ed2]DeU - e^[D„Sj{U)]ed,D,U 

+De (^-e^S2{U)[De,ed2]U - e^[Ds,Sj{U)]edjU^ . 

and in the same way, for 1 < i < 4, 

e'^dtiedifU + A{x, £Dx){edifU + e^Si{U){edif+^U + e^ S2{U)ed2{ediYlJ 

= -^^ Yl Ci{ediYSj{U)edj{ediY-''U . 

k=l 
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In all cases, the terms of higher order disappear by integration in x and the other terms are 
controlled with the following trilinear estimate (writing generically Q{U) for all the nonlin- 
earities): for aU < /c < 2 and all < £ < 4 

I < D';U\D^,Q{U) > I + I < {ediYu\{edifQ{U) > \ 

(7.5) ^ C\\U\\^i,^i\\D^ML^ + Il(e5i)^?7||z.2 + \\U\\L2f 

< j{\\dIu\\l2 + wiedi^uu^ + \\u\\L2y. 

Remark 7.1. Note that because of the bad embedding inequality \\\/U\\l°° < — ||C/||^^4, we 
lose one power of e, which seems not to be optimal considering for instance the fast oscillating 

j . A challenging question in order to apply semiclassical methods 

to nonlinear problems is to determine appropriate functional spaces (in the spririt of |23j ) 
which measures on the one hand the Sobolev regularity of the amplitudes, and on the other 
hand the oscillation frequency. 

We are finally able to obtain a uniform life span for the weakly nonlinear system, thus 
proving result (jlah of Theorem [2j Indeed combining the trilinear estimate (I7.5p and the 
commutator estimate (17. 3p . we obtain the following Gronwall inequality 



^'1 {WD'sUWh + \\i^d,)'U\\l, + \\U\\l,) < Ce' (l + WD'Ml^ + \\ie^^)'U\\L^. + \\U\\L^y 
from which we deduce the uniform a priori estimate 

WD^Mh + \\{edi)'U\\h + mh < {Co - Ct)-'' 

where Cq depends only on the initial data. Such an estimate shows that the life span of the 
solutions to (|7.ip is at least T* = Cq/C. 

7.3. Approximation by the linear dynamics. In this paragraph we shall prove re- 
sults (jlbp and ([2]) of Theorem[2j The proof of both results relies on standard energy estimates. 
We have 

e^dtiUe - Ve) + A{x,eD,){Ue - Ve) + e'^+'^Si{Ue)ediUe + e'^+'^S2iUe)ed2Ue = 
• If ?7 = and eV^ — t- in L°°, we use the decomposition 

e^dtiUe - Ve) + A{x, eD,){Ue - Ve) + e\Sj{Ue) - Sj{Ve))edjUe + 8'^Sj{Ve)8djUe = 

and obtain the following L^ estimate 

y^llf/. - K|li2 < 3e3||ea,?7,||Loo||f/, - Kill, +3e3||T4||^^||ea,-C/,||i2||?7, - K||i2 

< Ce\e\\edjUe\\L^ + \\edjUe\\l2)\\Ue -VeWl^ + Ce\e\\Ve\\L^)^ 
from which we conclude by Gronwall's lemma 

\\Ue-Ve\\l2<C f (e||K(s)||Loo)2expCr/ {eWedjUeh^ + \\ed,Ue\\l2)da] ds 

on [0, T*[, and that proves result (jlbp . 
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• If ry > 0, the same arguments show that the hfe span of the sohitions to (12. 5p tends to 
infinity as e — )■ 0: Tg > Ce"'' , and that these solutions are uniformly bounded in W^ on 
any finite time interval. Furthermore, on any finite time interval [0,T], the previous energy 
estimate gives 

y|r. - VeWh < Ce'+^ed,Ue\\L-\m\l, , 
from which we deduce 

\\U,-Ve\\l2<Ce^ [ e\\edjUeis)\\Lo^\\Ueis)\\l2ds . 
Jo 

Result ([2|) of Theorem [2] is proved. 

Appendix A. A diagonalization theorem 

In this appendix we shall state and prove the crucial theorem allowing to diagonalize semi- 
classically the matrix of pseudo-differential operators A{x2,£D). The construction of Rossby 
and Poincare modes (see Propositions 13.21 and 13.51 in Section [3]) are direct corollaries of that 
theorem. We refer to Appendix [B] for the notation and results of semi-classical analysis used 
in this paragraph. We consider an order function g on R^ , and a symbol h{X; 8, r) in S2d{g): 
depending polynomially on e and r, where we have defined X = {x,^) G R^ . The statement 
is the following. 

Theorem 3. Let g be an order function on 'R?'^ , and let h{X;e,T) be a classical symbol 
^^ S2d{g)! depending polynomially on e and t: there is an integer Nq such that 

No 

(A.l) VX G R^"', h{X-e,T) = Y^ h^'^{X)e^T\ 

j/=0 
where the symbols h^'^{X) belong to S2d{g)- Letfg = feiX) be a root of the polynomial h{X; e, •) 

oo 

which can be written for some z^ G N, fe(X) = > e'^'^ tj^{X) + 0{£°°), with tq 7^ and 

fc=o 
where tq is a symbol. Finally let /io(X, r) be the principal symbol of h, satisfying the following 

assumption: 

(A.2) 3C >0, yx e R^'^, Ve g]0, 1[, \drho{X; £''to)\ > C. 

Let K, be a compact subset ofH. Then there is a pseudo-differential operator T of principal 
symbol s^tq such that if ip € S'(R ) is microlocalized in K, and satisfies idgip = Tif) then 

(A.3) Op°(/i(e,i9,))V' = 0(e°°) in L^{Y{!^). 

Proof. The idea of the proof of the theorem is the following. Let us define a smooth function x, 
compactly supported in R^ , identically equal to one on /C. Then for any integer N , we shall 
compute recursively the coefficients of the symbol t^ := s'^tq+'^i^^-^ e'^^^t^ so T^ := Oy>\{t^) 
satisfies the required property, where r|^, unique up to 0(e°°), is given by 

00 

k=X 
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This allows to replace the root t^ by an actual symbol. The above strategy will be achieved 
in the following way. We notice that if idgip = T^ip then of course 

N 

idstp = Opl{T^'^)i; + 0(8^+"+^), where t^^^ := eVo + ^ e'^^'^XT^- 

fc=i 
If moreover ^ is microlocalized in /C, then by definition of x one has using (1B.3IIB.4]) and the 
fact that X is identically equal to one over /C idgip = Op^{T^ )iIj + 0{e '^^'^ ). This means 
one can (and shall) compute recursively {T^)i<k<N so that 

Opl{h{e,ids))ip = 0{e^+^), when z^.V = Op^(rf )V' + 0(e^+^). 

Note that it is convenient in the computations to compute T-^ as the "right" -quantization of 
the symbol t^. Now let us carry out the algebraic computations allowing to achieve the result. 
We shall start by dealing with the case when z/ 7^ 0, as the computations can be carried out in 
an easier way, and then we shall discuss the case when u = Q. Recalling that h{X] e, f£(X)) = 
we infer that h^'^ = if £ < i^ and 

(A.4) /iO'-(x,C) + /i^°(x,e)To(x,e) = 0. 

Now we recall that for any -0 in 5'(R ), one has 

The above integral, as all the ones appearing in this proof, is to be understood in the distri- 
butional sense. So if r^ and ip are such that 0'p^{Tf)il) = idgip with ip microlocalized in /C, 

(A.5) = (2vre)"'^/ e^^^^/i(x,e;e,Opi(rf))^(y)dyd^ + 0(e^+i). 

Now we need to compute h{x,£,;e,Opl{T^)). Using the fact that h is polynomial in r, 

h{x,^-e,Opl{T^)) = j;/i^-'V,e)H'(OpKrf))'' 

j/=0 
=: Opl{h{x,^,y,r];e)) 

where using compositions rules recalled in Appendix[Bl the symbol h{x, ^, y, rj; e) (in the {y, rj) 
variable) can be expanded as 

00 
h{x, e, y, r/; e) = e" J] hH^, C, V, ri)e^ + 0(e°°) 

k=0 

where the principal symbol of h{x^ ^, y, rj; e) is 

(A.6) e'^h^ix, e, y, r,) = e^ {h^^'^{x, + h'^^x, OMy, v)) ■ 

We recall indeed that according to Appendix IBJ one has [Op^{T^ )) =: Op^{'mj{y,ri;e)), 
where mi(y,r/;e) = r^ {y,r]) ^-i^d 

m,(y,7?;e)=J^^9,X-i(y>^;^)5^f(y,r/) + 0(e-). 

k>0 
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In particular the principal symbol of (Op^(re))"' is e'^^T^{y,ri), which yields (|A.6p . We notice 
that due to (|A.4h . (1A.6P imphes that 

(A.7) h^{x,^,x,O = 0, 

so that h{x,^,x,^]e) = 0{e'^^^). More generally, plugging the expansion of rrij into the for- 
mula defining h and noticing that h^''^{x, ^) = drhQ\T-=Q{x, ^), one finds that there is h^{x, ^, y, 77), 
depending on the symbol coefficients of h{x,£,), and in the (y,r/) variables on the sym- 
bols To{y, ri),T{ (y, r?), . . . , r^„i(y, rj) only, such that 

(A.8) VA;>1, h}'{x,i,y,ri) = h^{x,i,y,Ti) + drho\r=Q{x,i)T^{y,ri). 

Finally going back to (lA.Sh we find that 

Ov%h{E,\))^{x) = (27re)-2'^ / e''''^e''''^h{x,W,iT,e)i,{y')dydy'didr^ + 0{e''). 

We can first perform the integration in the y variable, which creates a Dirac mass at .^ — ry, 
and therefore we have 

(A.9) OpO(/i(e,A))V(x) = (2^e)-^ / e'^^'^h{x,i,y\i-e)i^{y')dy'di + 0{e''). 

Now we shall construct t^ ensuring that the order of Op^(/i(e, A))^(x) is 0{e'^~^^) instead 
of 0{e'^). The argument will easily be adaptable by induction to Op^(/i(e, A))V'(a;) = 0{e^) 
for any A^ G N, by a convenient choice of t~, for j < N — 1. We notice that in ()A.9p . 

the quantity h{x,^,y' ,^;e) can easily be replaced by h{x,^,x,^]e) by Taylor's formula: more 
precisely we write 

(A.IO) h{x, ^, y', C; e) = h{x, C, x, C; e) + [y' - x) • (Vy /i)(x, C, x, C; e) + 0{\y' - xf) 

which gives after integrations by parts 

OpO(/i(6,A))V(x) = / e'^''^h{x,C,x,C;e)^l^{y') '^^''^^ 



-is f V^e^'''^-{Vyrh){x,i,x,i;e)i^{y')^p^^+0{e-+^) 

A^-y')< ~ \ , , i\ dy'd^ 

e e h{x,i,x,e.;e)ip{y)— — r^ 

+ ie f e^'-^^V^.{v/h){x,i,x,i;e)i^{y')^^^+0{e''+^) 
Now using (IA.7p . it remains to choose t^ so that 

R2d •. « (27re)'' 

That is possible simply by looking at formula (jA.SP and choosing 



T"i {x,0 



dTho\r=o{x,0 
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Note that Assumption (|A.2|) guarantees that Tj" is weh defined. The argument may be pursued 
at the next order simply replacing (lA.lOp by 

+ ^{y'-x)^{y'-x):{Vl,h){x,^,x,^;e)+0{\y'-xf) 
and using integrations by parts again. Then the choice 

_ /i2(rE,e,x,e) + iV5-(Vj,/ii)(x,e,x,e)-iv2:(v2^o)(x,c,x,c) 

gives that Op^(/i(e, A))^(x) = 0{e'^~^'^). We leave the rest of the induction argument to the 
reader. To end the proof of Theorem [3] we need to consider the case u = Q. The argument is 
similar to the case z^ 7^ treated above, though the formulas are slightly more complicated. 
We shall use the same notation as in the previous case. We start by noticing that by definition 

of fe we have in particular for all (x,^) € R^*^, /J ^"'' (^^lO'^oC^'O = 0- Then as before let 

i=o 
us write h{x,^]e,0^i,{T^ )) =: Op^{h{x,^,y,r];e)). One computes easily that the principal 
symbol of h{x, ^, y, rj; e) is (unlike the case z^ 7^ 0) 

No 

(A.ii) h\x,ty,v) = J^h^'\x,Or^{y,v)- 

3=0 

Note that as above one has h^{x,^,x,^) = 0, so that h{x,^,x,^;e) = 0(e). One can also 
compute the next orders, and as in the case v ^ they can be written in the following form: 

Vfc> 1, h^ix,c,y,v) = h''ix,c,y,v) + rkiy,v)^h^'\x,0j4~\y^v) 

where h {x,^,y,r]) depends on the symbol coefficients of h{x,^), and in the (y,r/) variables 
on ro(y, rj), t^ {y, 77), ... , TJ^_-^^{y, rf) only. In particular we notice that 

(A.12) VA:>1, P(x,e,x,e) = /i'=(x,^,x,e)+r-(y,r7)5,/io|.=.o. 

Now that these formulas have been established, it remains to go through exactly the same 
computations as in the case i^ = 0, and we find 

_ h^{x,i,x,i)+i'SIs,-{^y^G){x,i,x,i) 

which is well defined thanks to Assumption (|A.2p . The other orders are obtained exactly as 
in the case z/ = 0. This ends the proof of Theorem [3l D 

Appendix B. Some well-known facts in semi-classical analysis 

In this section we recollect some well-known facts in semi-classical analysis, which have been 
used throughout the paper. Most of the material is taken from [7], |21] . [25], |30] and [29]. 

B.l. Semi-classical symbols and operators. 
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B.1.1. Definitions. We recall that an order function is any function g G C°°(R ; R^ \ {0}) 
such that there is a constant C satisfying 

yX G R'^, Va G N^ |9°5(X)| < Cg{X). 

For instance g{x,^) = (l + |i^p)2 =: {£,) is an order function. Note that the variable X usually 
refers to a point {x, ^) in the cotangent space r*R" = R^", or to a point of the type {x, y, ^) 
with y G R". A semi-classical symbol in the class Sfi{g) is then a function a = a{X;£) 
defined on R x]0,eo] for some Eq > 0, which depends smoothly on X and such that for 
any a G N*^, there is a constant C such that \d°'a{X,e)\ < Cg{X) for any {X,e) G R'^x]0, eo]- 

If (aj)jgN is a family of semi-classical symbols in the class Sdig), we write that 

oo 

i=o 
if for any A^ G N and for any a G N'^, there are eg and C such that 

N 

yX G R"', Ve G]0,eo] d''(a{X,e)-^e^aj{X,e)] < Ce^g{X). 

j=0 

Conversely for any sequence (aj)jgN of symbols in Sd{g), there is a G Sd{g) (unique up 

oo 

to 0(e°°)) such that a = V^e-'aj + 0(e°°). An e-pseudodifferential operator is defined as 

i=o 
follows: if a belongs to Ssnio), and u is in P(R'^), then 

Op,(a))n(x) := ^-L^ J e^(-~y)</-aix,y,Ou{y) dyd^. 

B.1.2. Changes of quantization. Ifa G S2n{g) and t G [0,1] then a*(x,y,^) := a{{l—t)x+ty,£) 
belongs to S^nid), and one defines Op* (a) := Opj(a*). When t = this corresponds to the 

classical, or "left" quantization, when t = 1/2 this is known as the Weyl quantization 

1 

(and is usually denoted by Op^(a) = Op|(a)), while when t = 1 one refers to the "right" 
quantization. Furthermore, if a belongs to S^niiO'^) fo'^ some integer m (or more generally 
if a belongs to S^nid) where 5 is a Hormander metric [2l|), then there is a unique symbol at 
belonging to 5'2n((0™) (resp. at G Ssn{g)) such that Op*(at) = Op^{a), and one has 



at(x,e) = j; ^^I'S.Ma^ + t7?,x - (1 - t)7?,e)|,=o + 0(^°°)- 



A classical symbol is a symbol a in S2n{{0"^) such that a{x,^;e) = y^e-'oj(x,^) + 0(e°°) 

i=o 
with flQ not identically zero, and a-,- G 5'2n((0"^) independent of e. For any real number v, the 
term e'^a^ is the principal symbol of the classical pseudo-differential operator A = e'^Op*(a) 
(and this does not depend on the quantization). On the other hand e'^'^^ai is the subprincipal 
symbol of ^ = e''Op^(a) (in the Weyl quantization only). In the following we shall denote 
by a't{A) the symbol of an operator A = Op* (a) (in other words a = at{A)), and by ap{A) 
its principal symbol. 
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B.1.3. Microlocal support and e-oscillation. Ifu is an e-dependent function in a ball of L^(R"), 
its e-frequency set (or microlocal support) is the complement in R^" of the points (xo,i^o) 
such that there is a function xo £ 5'2n(l) equal to one at (xc'^o)) satisfying 

l|Opf(xo^)||L2(R2n) = 0(e°°). 

We say that an e-dependent function /^ bounded in L^(R") is e-oscillatory if for every 
continuous, compactly supported function (p on R"", 

(B.l) limsup / \ipfe{Ofd^-^0 as R-^oo. 

An e-dependent function /g bounded in L^(R'^) is said to be compact at infinity if 
(B.2) limsup/ \fe{x)fdx^O as R ^ co. 

e-5>0 Jx>R 

B.1.4. Adjoint and composition. Let a be a symbol in S^nid), where g is a Hormander met- 
ric [21], and define a*{x,y,£,) := a{y,x,^). Then the operator {Op^{a))* := Op^{a*) satisfies 
for all u,v in 5(R"'), 

1 
and is therefore called the formal adjoint of Opi,{a). In particular Op| (a) is formally self- 
adjoint if a is real. Let a and b be two symbols in S2n{gi) and S2n{92) respectively, where Qj 
are Hormander metrics. For all t € [0, 1], there is a unique symbol q in S2n{9i92) which allows 
to obtain Op* (a) o Op* (6) = Op*(q). Moreover one has 

(B.3) ct{x, e; e) = g-I^-^e-^.^.] (^((^ _ ^^^ ^ ^^^ ^)^((i _ ^^^ + tx,0) iu = v = . =■ ^#'b. 

This can be also written 

ctix, C; e) = E ^(^"^^ - ^?^")' (°((^ - *)^ + *^' ^)^((^ - t> + te, 0) I , = „ = , + 0(e°°). 

fc>0 ■ r) = ^ 

In particular one has at{A o B) = at{A)at{B) + 0{e). For example in the case when i = 
then Op£(a) o Op^(b) = Op^{c), with 



(B.4) c{x,0 = a#b=^^dfa{x,0d^bix,0 + Oie^). 



In particular if a and b are two classical symbols in the sense described above, and if one 
defines A := Op£(a) and B = Op^{b), then the principal symbols satisfy (if ap{A)ap{B does 
not vanish identically) ap{AB) = ap{A)ap{B). 

B.2. Semiclassical operators, Wigner transforms and propagation of energy. One 

of the main interests of the semiclassical setting is that it allows a precise description of the 
propagation of the energy, on times of the order of 0(e). We refer for instance to [H] (Section 
6) for the proof of the following property (actually in the more general setting of matrix- valued 
operators): consider a scalar symbol Ti;{x,^) defined on R^", belonging to the class S2n{{0'^) 
for some cr G R (or more generally to S2n{g) where g is a Hormander metric). We assume 
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moreover that Op^(r£) is essentially skew-self-adjoint on L^(R"). Then consider /^ an e- 
oscillatory initial data in the sense of (IB.ip . bounded in L^(R"') and compact at infinity in 
the sense of (JR2i) . and the PDE 

edtfe + Ovf{re)fe = ^. fe\t=0 = fe- 

Then the Wigner transform W!r{t,x,(,) of fe{t) defined by 

converges, locally uniformly in t, to the solution W of dtW + {to,W} = where tq is the 
principal symbol of r^, and where the Poisson bracket is given by 

{to, W} := V^To • V^W - V,.To • V^W. 

The interest of Wigner transforms lies in particular in the fact that under the assumptions 
made on /^, for any compact set K C R" one has J^ [/^(t, x)p dx = W£{t,K x R") due to 
the fact that \fe{t,x)\'^ = Jj^„ We{t,x,^) d£,. 



B.3. Coherent states. A coherent state is ^p^q{y) := (7re)~4e* J e 2^ . Any tem- 
pered distribution u defined on R" may be written 



u{y) = (27re) 2 / Tu{p,q)<^p,q{y) dpdq, 
where T is the FBI (for Fourier-Bros-Iagolnitzer) transform 

n 3n / ■(l — y)-p (y — Q) 

Tu{p,q) :=2~2(7re) 4 / g* J e ^^u{y)dy. 



This transformation maps isometrically L^(R") to L^(R^"). The above formula simply trans- 
lates the fact that u = T*Tu. 



B.4. Fourier Integral Operators. A Fourier Integral Operator (FIO) is an operator U 
which can be written, for any / G L^(R") 

Uf{x) = —^ f e**(^'^'^)/^a(rE, y, T)f{y) dydr 
(27re)~ JR2" 

where a is a symbol of order 0, compactly supported in x and y, $ is real valued and homoge- 
neous of degree 1 in r, smooth for r ^ 0. One requires also a non degeneracy condition on the 
phase (see [29], Chap. 9, Par. 6.11) on the support of a. Then U is continuous over L^(R'^). 
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